Non-perturbative renormalization group for the Kardar-Parisi-Zhang equation 
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We present a simple approximation of the non-perturbative renormalization group designed for 
the Kardar-Parisi-Zhang equation and show that it yields the correct phase diagram, including the 
strong-coupling phase with reasonable scaling exponent values in physical dimensions. We find 
indications of a possible qualitative change of behavior around d — 4. We discuss how our approach 
can be systematically improved. 
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The Kardar-Parisi- Zhang (KPZ) equation Q: 
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dh(x, t) 
dt 



v\J 2 h{x,t) + - (Vh(x,t)) 2 + rj(x,t) (1) 



where rj is an uncorrelated Gaussian noise with zero mean 
and variance (r](x, t)r](x',t')) = 2 DS d (x — x') S(t — t') is 
maybe the simplest nonlinear Langevin equation show- 
ing non-trivial behavior 0]. Under form (p} it describes 
the kinetic roughening of a d-dimensional interface [ij, 
but it is now recognized as representing an extremely 
large class of non-equilibrium or disordered systems [2j . 
A large body of works at all levels have led to a general 
agreement about the phase diagram of the KPZ equa- 
tion, usually described in terms of the dynamical scaling 
properties of h(x,t), characterized by the two-point cor- 
relation function C{\x-x'\, t-t') = ([h(x,t)-h(x',t')] 2 ). 
In particular, at large scales, C shows the scaling form 
C(L,t) = L? x J(t/L z ) with the two exponents related 
via z + x = 2. For d < 2, the nonlinear term is rel- 
evant and the interface always roughens (x > 0). For 
2 < d < 4, two phases exist depending on the value 
of the effective coupling constant g = \ 2 D/v 3 : in the 
strong-coupling regime, rough interfaces are again ob- 
served, while for small g values the nonlinear term is 
irrelevant and the interface is smooth (x = 0, z = 2). 

Some important points remain controversial: For d > 
4, numerics Q and real-space calculations Q still show 
the existence of both a rough and a smooth phase, while 
some theoretical approaches 0-0] argue that d = 4 is 
some kind of upper critical dimension beyond which 
2 = 2. (For a recent discussion, see Q.) Thanks to an in- 
cidental extra symmetry, x — h is known in d = 1 Q , but 
no exact results are available in higher dimensions. This 
unsatisfactory situation is largely due to the fact that the 
strong coupling phase of the KPZ equation has remained 
out of reach of controlled analytical approaches. In par- 
ticular, standard perturbation expansions fail at all order 
to find a strong coupling fixed point [ltj. Some non- 
perturbative approaches have been deviced, such as the 
mode cou pling approximation 0, [Tljj , the self-consistent 
expansion [l2j . or the weak noise scheme 0] , but they are 



difficult to improve in a systematic way. In this context, 
the non-perturbative renormalization group (NPRG) [l3| 
appears promising, in particular since it is able to deal 
with perturbatively unaccessible fixed points even in out- 
of-equilibrium problems 14]. In this Letter, we present 
a simple approximation of the NPRG designed for the 
KPZ equation and show that it yields the correct phase 
diagram, including the strong-coupling phase with rea- 
sonable scaling exponent values in physical dimensions. 
In particular the strong-coupling fixed point in d > 2 is 
genuinely non-perturbative (hence essentially out of reach 
of perturbative RG) and fully attractive — which roots 
the existence of generic scaling. We find indications that 
a qualitative change of behavior occurs around d = 4, but 
at the minimal order presented here, our approximation 
does not resolve the puzzle of the existence of an upper 
critical dimension. We finally discuss how our approach 
can be systematically improved. 

We first recall how the KPZ problem can be cast into 
a field theory 15]. Introducing a response field h and 



sources the generating functional reads: 

2[j,j] = J V[h, ih] exp f-S[h, h] + J(j h + jh)\ (2) 

S[h,h] =J Ui (d t h-vV 2 h-^{\Jh) 2 ^j -D/?j(3) 

where x = (x,t). The symmetries of the KPZ equa- 
tion are three-fold: (i) invariance of Z under the 
Galilean transformation 7g = {M x ) ~ * + ^vt,i) + 
v.x,h(x) —¥ h(x + Xvt,t)}] (ii) invariance of the com- 
bination S — J hdth under the "time-gauged" symme- 
try Ttg = {ft(x) /i(x) + f(t)} where f(t) is an 
arbitrary function of time; (iii) in = 1, additional 
time-reversal invariance 7tr = {h(t) — > —h(—t),h(t) — > 
h(-t) + ^V 2 h(-t)}. 

The NPRG builds a one-parameter family of models 
indexed by a scale k such that fluctuations are smoothly 
included as k is lowered from the microscopic scale A, 
where no fluctuations are taken into account, to k = 
where they have all been summed over [13] . To this aim, 
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we add, to the original action S, a momentum-dependent 
mass-like term AS k which preserves the symmetries: 
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AS k = 7: I hi [Rk]ij hj ; Rk=r{-^ 





Vkq 2 



-2D k 



(4) 



where q = (q,uj), q = \\q\\, i,j <G {1,2}, hi = h,h 2 = h, 
and summation over repeated indices is implicit. With 
the choice r(x) = a/ (exp(x) — 1) where a is a parameter, 
the fluctuation modes hi(q > k) are unaffected by ASk, 
while the others (hi(q < k)) are essentially frozen. The 
effective action V'Jj where ipi — (hi) are the expec- 

tation values of the fields hi in presence of the external 
sources j, j, is given by the Legendre transform of log Zk 
(up to the term proportional to R k ) 13]: 

T k [ip, ip] + log Z k [j, j] = [jiipi ~ \ I i)i [Rk]ij ipj (5) 



From T k , one can derive 2-point correlation functions 

s 2 r k [i>,ip] 



-,(2)- 



! (xi,x 2 ,V,V') 



(6) 



Sip tl (yii)S^i 2 (x 2 ) 

and more generally n-point correlation functions that we 
write, for future convenience, in a 2 x 2 matrix form as 
(omitting the dependence on the running scale fc) 

^-2 r (2)( Xl)X2i ^) 



-.(") 



His ( X 3)---^i„( x n) 



(7) 



The exact flow for T k [ip, ip] is given by Wetterich's equa- 
tion El: 



fc r fc = iTr/<v, /?, -G,, with G,= 



,(2) 



Ri 



(8) 



When A: flows from A to zero, T k interpolates between 
Tfe = A = S and the full effective action T k =o- Solving 
Eq.© is thus equivalent to solving the model. This is 
however impossible to do exactly because it is a non- 
linear integral partial differential functional equation. 

Here, as usual, our interest lies in the fixed point struc- 
ture, which is determined by the small momentum sec- 
tor. However, the signature nonlinear term of the KPZ 
problem involves gradients, and this seems to re quir e the 
faithful description of the full momentum range [16| • In- 
deed, expanding the rj^'s in powers of their momenta 
(the so-called derivative expansion), yields uncontrolled 
results at least at low orders [l6[. Here, by contrast, we 
define a strategy which yields closed flow equations for 
the correlation functions, and detail it for the two-point 
functions r^ 2 '. 

We first differentiate Eq.© twice with respect to the 
fields and evaluate it, without prejudice for the following, 
in uniform field configurations ^ u = {?/>, V'lunif : 



fc [r< 2 >]y(p) = Tr^9 fc i?(q).G(q)^-irW(p,- P! q) 



+rf ) (p,q).G(p + q) 



rf(-P,P- 



q) • G(q) (9) 



where the k and dependences have been omitted, as 
well as the last argument of the T^ n \ thanks to trans- 
lation invariance. Our approximation closes the exact 
Eq. by approximating the three- and four-point func- 
tions. To this aim, we make use of the following exact 
relations, valid in uniform field configurations: 

4 m '")(q 1 =0 ) {q i } ) *u) = ^ri m - 1 '")({q i },* u ) (10) 
ri m < n) ({qa,q m+ „=0,* u ) = ^rf^lqj,^,) (11) 

where, e.g., I^ 2 ' 1 ) is understood as [T]?'] ,7, and more 
generally r' m '") is the r( m+n ) function involving m (resp. 
n) derivatives w.r.t. tp (resp. "0). To set one q to zero, we 
use the following facts: Thanks to the presence of the cut- 
off function R k , the momentum dependence of the corre- 
lation functions is smooth and the internal momentum q 
in ((9]) is effectively limited to q = \\q\\ < k. Setting q = 
for the rj? and T^ 4 '' functions in © all along the flow 
is thus a reasonable approximation which becomes exact 
for any finite external momentum in the k — > limit. An 
exception, though, must be made for the F^ 2,1 ) function: 
its bare momentum dependence Xq\ ■ q\ — which stands 
for the signature nonlinear term of the KPZ problem — 
must be kept and the approximation above only imple- 
mented on the other terms of the function. In addition, 
we restrict ourselves to zero external frequencies. Finally, 
to be able to use (flQl) and (fTTj) , we set the internal fre- 

(see below). At 



quencies to zero in rjj, 3 ' and in 



this point, the and T k ^ are replaced by derivatives 
of the T k and Eq. (j9|) is closed. In fact, the symmetries 
impose that d^,T k n ' = 0, because of (fl"0|) and of the Ward 

identity T k m ' n \q x = 0, Wi,q 2 , . . . ,q m+ „-i) = (derived 
from 7tg and valid for m > and (m, n) 7^ (1, 1) [171]). 

Solving the problem at this level (with one field de- 
pendence) was done in equilibrium theories with great 
success [18|. Here, this remains very difficult because the 
KPZ symmetries are not easily maintained along the RG 
flow. Thus, on top of the approximation above, we also 
perform an expansion around ip = and keep only the 
leading terms present at the bare level. All T k and T^ 4 - 1 
then vanish, except T^ 2 ' 1 ) = Xqi ■ q\. Using various Ward 
identities imposed by the KPZ symmetries, we find that 
the two-point functions are completely determined. In 
particular, 7g and 7tg impose that the bare frequency 
dependence is not renormalized, so that, finally, rl 1 ' 1 ^ 
and rj. ' 2 ^ are parameterized by two arbitrary functions 
7fc 1 (<?) and 7° 2 (g), yielding our Ansatz [171 ]: 



r 



-.(1,1) 
-.(0,2) 

k 

[2,0) 



(q,uj,tp) = iuj + q 2 jl 1 (q) 

(q,u,f) = iFill 
(q,u,ip) = -\q 2 4> . 



(12) 



Note that, in d = 1, where the extra time-reversal sym- 
metry holds, 7° 2 = —{2Dk/vk)l\ 1 - Coming back to the 
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internal frequency u, our Ansatz allows to compute ex- 
actly the integrals over uj (which we find convergent). 
Thus, the large-w domain does not contribute much, 
which justifies a posteriori our neglecting of higher-order 
terms in ui. In fact, we also set to zero the term iuiq 2 in 
■^k > w hich is small for the momenta of interest. 

To treat efficiently the zero-momentum sector of in- 
terest near a fixed point, we introduce dimensionless 
and rcnormalized quantities. All momenta are mea- 
sured in units of k (e.g. q — q/k). At the bare level, 
7 A 1 (<7) = v and 7 A 2 (g) = —2D. We thus define the di- 
mensionless renormalized functions ji (q) = 7fc 1 (<z)/^fc 
and 7° 2 (<z) = The running anomalous di- 

mensions T] v (resp. rj D ) are defined by kdkVk = —flu{k)vk 
(resp. kdkDk = —r) D (k)Dk) so that at a fixed point, 
Vk ~ k~ Vv and Dk ~ k~ VD . The scaling exponents are 
then expressed as z = 2 - 77* and \ — (2 - d + rf D - rj*)/2. 
To fix the absolute normalizations of 77: (q) and (resp. 
7° 2 (g) and D k ), we set 7^(0) = 1 (resp. 7 ° 2 (0) = 1). 
(This choice is dictated by the fact that the Ansatz is 
designed around the small-q sector.) 

Inserting the Ansatz into ©, we checked that A does 
not flow (kdkF^' ^ = 0), in agreement with the Galilean 
invariance, and thus the dimensionless running coupling 
constant reads g k = fc d ~ 2 A 2 Z?fc/i/|. Its flow equation is: 



d s gk = 9k(d - 2 + 3r)„(k) - i] D (k)) , 



(13) 



where d s = kdk, so that z + \ = 2 is enforced at any 
fixed point with g* ^ 0. We are finally left with the 
following flow equations for the two running functions 
(omitting the k index, using u = p+q, P = —p 2 +q 2 +u 2 , 
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FIG. 1: 
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(Color online) (a) Flow of the 7 (p) function from 
its bare shape (7° 2 (p) = 1) to the KPZ strong coupling fixed 
point Fsc (solid red line) in d = 3. (b) Flow of gvd from 
various initial bare values in d — 3. (The normalization con- 
stant v' 1 = 2 d+1 7r d/2 r[f] is related to the integration vol- 
ume). Flow lines in red (resp. purple) converge to Fsc (resp. 
-Few). The unstable fixed point Frt is given by the blue 
dashed line, (c) Flow diagram within our approximation in 
the (d, gvd) plane. Red circles: renormalized value gg C at 
Fsc- Dashed purple line: Gaussian Few fixed point. Cyan 
squares: bare value g c separating the basins of attraction of 
Fsc and Few- Grey lines symbolize flow lines, (d) Variation 
with d of x = 2 — z for Fsc (red circles: our results; orange 
squares: numerical values from [1,3- See also Table UJ) 



d s f\p) = {ilo+pd p )^{p)-~g -= 



9 s 7 11 (p) = (^+pap)7 11 (p)+.9 



P i ~g 02 (u)H(q,u) 



L(q,u) 



p2 



L(q,u) 



(14) 



U 2 - Q 2 \ 

-=2 9ii{u)H(q, u) - —g u (q) g 02 (u) r v (q) j (15) 



TABLE I: Exponent values in integer dimensions. The aver- 
age numerical values are extracted from 0]. To our knowl- 
no estimates of uj are available in the literature. 



d 


1 


2 


3 


4 


X (this work) 


0.50 


0.33 


0.17 


0.075 


X (numerics) 


0.50 


0.38 


0.30 


0.24 


uj (this work) 


0.817 


0.70 


0.63 


0.54 



with H(q,u) = gii(q)((gii(q) + gii{u))r D {q) 

-502(g) (2511(g) +fci{u))r v (q) (16) 

and L{q,u) = 5ii(g)5ii(«)(.9ii(<?) + .9n(«)) 2 (17) 

where r v (q) = rj v r(q) + qr'(q), gn{q) = g 2 (7 n (<7) + r(q)) 
and similar definitions for td and 502- (Note that in 
d = 1, the 7 02 and 7 11 functions are identical.) 

We have performed the numerical integration of 
Eqs. (I13|14|15|) (together with the renormalization condi- 
tions 7 02 (0) = 7 11 (0) = 1) by discretizing momentum on 
a ^fp mesh of typical spacing 5.10~ 2 , using 5-point finite- 
difference expressions for the pdp terms, Euler explicit 



time-stepping with a typical time step As = — 10~ 4 . All 
integrals were estimated using Simpson's rule and Carte- 
sian coordinates. We have checked the robustness of the 
results presented below against numerical resolution and 
against variations of the momentum range used for the 
cut-off function r (typically 4 in units of p). In short, 
the simplest numerical techniques performed nicely, with 
typical runs taking minutes on a current computer. 

We systematically followed the flow of functions 7 02 
and 7 11 , together with the running coupling constant g 
and running exponents r\r> and r\ v from the bare initial 
condition 7° 2 (p) = J 11 (p) = 1 down to k = (s — > —00). 
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Figure QJ summarizes the flow diagram obtained in our 
approximation. In all dimensions studied (i.e. up to 
d = 8), we have found, besides the Edwards- Wilkinson 
fixed point i^Wi a fully-attractive non-trivial fixed point 
Fgc- In ah d, we find generic scaling, i.e. the flow al- 
ways reaches one of these fixed points. For d < 2, Fgc 
is reached from any initial condition. Along the flow, 
the 7 02 (p) and 7 (p) functions gently deform to reach a 
fixed shape with an algebraic tail at large p (Fig.[T^). For 
d > 2, Few becomes locally fully attractive, and there 
exists a critical bare value g c separating the basins of at- 
traction of Fgc and Few (Fig. QJ>). Right at g c and for 
d < 5, the flow reaches an unstable fixed point Frt which 
drives the roughening transition. In d = 2, Frt coincides 
with Few, and becomes non-Gaussian for larger dimen- 
sions. Because Galilean symmetry is respected (and its 
coupling constant <?g C is non zero), Fgc is characterized 
by a single exponent. We thus only discuss below the 
values obtained for x 0- Our work also gives access to 
the subleading exponent uj characterizing the approach 
to Fgc- Table U and Figure QJl contain our estimates. 
For d < 4, x decreases almost linearly with d, with the 
exact value \ (resp. §) recovered in d — I (resp. d = 0), 
and a reasonable but deteriorating agreement with nu- 
merical values for 2 < d < 4. In higher dimensions, 
X increases with d, at odds with both numerical values 
and the scenario of d = 4 being an upper critical di- 
mension beyond which x = [5H7|j. Regarding Frt, we 
record negative values of x f° r 2 < d < 5, which is rem- 
iniscent of perturbative results performed at fixed d [2(J 
but in contradiction with exact results dictating x = 
Tol 21 1 . For d > 5, Frt seems unreachable from bare 
initial conditions, probably an effect of the crudeness of 
our approximation. 

Several observations hint at a qualitative change occur- 
ing around d — 4: the values of gg C and g c first increase, 
then decrease (Fig. [Tfc). For large d, the unstable fixed 
point Frt becomes ill-defined. At this stage, we are un- 
able to decide whether the change of behavior near d = 4 
is related to the existence of an upper critical dimension 
or whether this is just an effect of the deterioration of 
the quality of our scheme in large d. Higher orders of our 
approximation, if only by adding the next terms in ip to 
the two-point functions, might yield the answer to this 
important question. At any rate, we believe our results 
represent a breakthrough, essentially because our approx- 
imation scheme is able to deal with correlation functions 
in a non-perturbative and systematically improvable way. 
In particular, as opposed to the mode-coupling approx- 
imation, which assumes scaling forms for the two-point 
functions, our equations flow generically to a fully-stable 
strong-coupling fixed point — which roots the existence 
of generic scaling — , or to the EW fixed point for d > 2, 
yielding, for the first time, the correct phase diagram of 
the KPZ equation within a RG approach. 
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